Self-consistent field theory for the interactions between keratin intermediate filaments 
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We have developed a model for the interactions between keratin intermediate filaments based 
on self-consistent field theory. The intermediate filaments are represented by charged surfaces, 
and the disordered terminal domains of the keratins are represented by charged heteropolymers 
grafted to these surfaces. We estimate the system is close to a charge compensation point where 
the heteropolymer grafting density is matched to the surface charge density. Using a protein model 
with amino acid resolution for the terminal domains, we find that the terminal chains can mediate 
a weak attraction between the keratin surfaces. The origin of the attraction is a combination 
of bridging and electrostatics. The attraction disappears when the system moves away from the 
charge compensation point, or when excess small ions and/or NMF-representing free amino acids 
are added. These results are in concordance with experimental observations, and support the idea 
that the interaction between keratin filaments, and ultimately in part the elastic properties of the 
keratin-containing tissue, is controlled by a combination of the physico-chemical properties of the 
disordered terminal domains and the composition of the medium in the inter-filament region. 

PACS numbers: 87.15.A-, 68.47.Pe, 05.70.Np 



I. INTRODUCTION 

The outermost layer of skin, the stratum corneum 
(SC), is often described as organised into a 'bricks-and- 
mortar' type structure, where the mortar represents the 
self-assembled lipid lamellae and the bricks refer to the 
protein-rich corneocytes. [iHl] • Corneocytes are nonvi- 
able disk-shaped flat horny cells mainly composed of 
keratin proteins, organised in complex intermediate fil- 
ament (IF) networks. Keratins, in turn, are important 
structural proteins which confer stiffness to many biolog- 
ical tissues such as skin, nails and hair. There are 54 
functional human keratin genes, of which 28 are type I 
(acidic) keratin genes and 26 are type II (neutral and ba- 
sic) keratin genes. A new systematic nomenclature and 
functional role of keratins was presented by Schweizer et 
al. [5], Moll et al. [Q , and Gu and Coulombe [7 . 

Keratin monomers consist of central a-helical rod do- 
mains of similar substructure (sa 310 amino acids) and 
two disordered (unstructured) glycine-rich N- and C- 
terminal domains of variable size. Two keratin polypep- 
tides associate in a parallel arrangement to form an 
« 50 nm long coiled coil dimer, consisting of two differ- 
ent type of keratins: one acidic (type I) and one neutral- 
basic (type II). The most frequent keratin (K) dimer ex- 
pressed in the SC and the upper epidermis is the Kl/KlO 
pair [5HTU] . The two coiled-coil heterodimers further self- 
assemble into tetramers by packing into an antiparal- 
lel staggered configuration. Tetramers, in turn, aggre- 
gate end-to-end forming protofilaments with a diameter 
around 2-3 nm. Two protofilaments make a protofibril 
with diameter of order 4-5 nm; four of these assemble 
laterally to form the keratin IF with diameter of order 
8-10 nm [3 [TOHTS] . Schematically, a keratin IF could 
be pictured as a long cylindrical object filled mainly by 



a-helical coiled coils domains, and decorated on the sur- 
face by disordered N- and C-terminal domains extending 
into the surrounding solution. An illustration of the IF 
hierarchical organisation is depicted in Fig. [T] 

Inside the corneocytes IFs are surrounded by a com- 
plex mixture of water, ions, free amino acids and other 
low molecular weight water soluble non-ionic compounds; 
this mixture is sometimes referred to as the "Natu- 
ral Moisturising Factor" (NMF). NMF plays an impor- 
tant role in skin moisturisation and in maintaining the 
physico-chemical properties of the skin, such as elastic- 
ity and permeability |19H261; and it results from pro- 
teolytic degradation of filaggrin, a histidine-rich protein 
[TOH^ . A reduced amount of NMF correlates with dry, 
flaky and itchy skin. Dry skin conditions may be a cos- 
metic problem triggered by natural (seasonal) changes of 
SC physical properties [501 HI]; but they may escalate 
to severe inflammatory skin disorders such as atopic der- 
matitis [22 [21], xerosis [22 [53] ichthyosis [H ^ and 
psoriasis [24l [26] . 

Jokura et al. [T^ studied the effect of NMF on SC 
elasticity using NMR spectroscopy, rheology and elec- 
tron microscopy. The authors observed that treating an 
excised SC sample with water releases NMF and leads to 
a keratin IF mobility reduction, and overall corneocyte 
rigidity. Electron micrograph evidence suggested that in 
the absence of NMF, keratin filaments tend to associated 
more tightly with each other. Further hydration of the 
sample does not improve the mobility of the fibers. How- 
ever, the original IF mobility conditions were partially 
restored by application of amino acid solutions. The au- 
thors compared the effect of different types of amino acids 
on the restoration of the SC elastic properties: neutral 
or basic amino acids, such as glycine or lysine, provided 
remarkable recovery of SC elasticity. In contrast, acidic 
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FIG. 1: IF organisation and the SCF model of N- and C- 
terminal domains attached onto IF surfaces. The IF surfaces 
are modelled as plane walls, the grafted domains as connected 
monomers, the salt ions and/or free amino acids as single 
monomers. All the other space is occupied by water. The 
separation between walls is varied in order to obtain the inter- 
action potential mediated by the walls with grafted domains. 



amino acid, such as aspartic acid, was not as effective. 

These findings suggested the hypothesis that loss of 
SC elasticity is due to increased intermolecular attractive 
forces between keratin filaments. In physiological condi- 
tions, NMF plays the important role to reduce these at- 
tractive forces, and to ensure SC elasticity. It is tempting 
to argue the protruding non-helical regions (unstructured 
N- and C- domains) mediate the interaction between the 
NMF-rich matrix and the IFs. In this work we present a 
modelling study of the interactions between keratin IFs 
suspended in different media: (i) a salt free solution mim- 
icking the NMF depleted system, and the effects of (ii) 
added salt and (iii) NMF-rich amino acid solution. 



II. METHOD 



A. Self-consistent field (SCF) approach 



FIG. 2: Schematic illustration of the N and C unstructured 
terminal domains for the AA and PG models. 



Interactions between two IFs formed by Kl/KlO ker- 
atins were investigated using the lattice self-consistent 
field (SCF) method [27 33 . The hehcal cores of the two 
IF are modelled as planar surfaces at distance D apart 
with disordered N and C terminal domains uniformly 
grafted onto them. The space between the IF surfaces is 
filled by water molecules, ions and/or free amino acids. 
The schematic model system is illustrated in Fig. [l] 



In the lattice SCF scheme the space between two sur- 
faces is divided into layers z — 1,2,3, D parallel to 
the walls, and each layer is further divided into lattice 
cells of equal size. Each lattice site is occupied by one 
of the monomeric species of the system (i.e. by pro- 
tein residue, water molecule, ion, etc.), so the total vol- 
ume fraction for all the species in each layer equals one. 
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= 1, where the volume fractions 4'°'{z) have 
the meaning of dimensionless concentration of species 
type a at the distance z from the surface. Obtaining 
the equihbrium concentration profiles for all the system 
components, is the primary target of the SCF 

calculations. The volume fraction distributions depend 
nonlinearly on the potential of mean force, act- 
ing on each species a in the system. The potential for 
each component, ^"(z), in turn, depends on the volume 
fraction profiles, as well as on the short range (Flory- 
Huggins) and long-range (electrostatic) interactions be- 
tween all the species of the system. To find both quan- 
tities, (/>"(2) and u"{z), a set of nonlinear equations is 
constructed and solved self consistently by an iterative 
procedure. The volume fraction profiles obtained in this 
way minimise the free energy of the system |31j . 

The SCF method is widely used to study properties 
of disordered proteins at interfaces. Earlier, the scheme 
was implemented to investigate adsorption of milk pro- 
teins, /3-casein [31] and asi-casein [55H55] . More recently 
the SCF approach was applied to protruding terminal do- 
mains of neurofilaments (NF) j39 -43j and to microtubule- 
associated 3RS tau protein, expressed in neurons of the 
central nervous system The detailed description of 
the method can be found in the original literature. Here 
we apply the method to study unstructured terminal do- 
mains of skin keratin IF. We have considered and com- 
pared two models for terminal N and C domains, detailed 
next. 



B. Amino acid (AA) model for terminal domains 

The first model (AA) is based on the primary struc- 
ture of N and C terminal domains for the keratins Kl 
and KIO. The amino acid sequence of these domains is 
taken from the Human IF Database |1S]. The terminal 
domains of Kl (Nl and CI) consist of 180 and 151 amino 
acids, respectively, and of KIO (NIO and CIO) of 146 
and 124 amino acids. All the amino acids in this model 
are divided into five groups according to their properties: 
'H', hydrophobic (Ala, Val, Leu, He, Met, Trp, Phe, Pro, 
Cys); 'P', polar (Ser, Thr, Tyr, Asn, Gin), 'G' (Gly); 
'-|-', basic (Arg, Lys, His); and '— ', acidic (Glu, Asp). A 
similar approach for allocating amino acids into groups 
is widely used in literature [34U44j . 

As can be seen from the amino acid sequence, N and C 
domains for Kl and KIO are glycine-rich (~ 50%), where 
glycine is mostly accumulated in blocks of 3-6 residues 
separated by one or two H or P residues. Being aware 
that glycine is a peculiar amino acid, showing both polar 
and hydrophobic behaviour (depending on the length of 
the poly-glycine residue) [15H49j . we reserve for glycine 
a separate classification group 'G'. 

The value of pH for SC is varied depending on SC 
depth, location and environment and the reported values 
of SC pR are in the order of 5-7 According to 

pKq values, at pH = 7 the amino acid residues Arg and 



Lys have charge q = +1 e. His has q = +0.36 e, and Glu 
and Asp have q — —le. At pH = 5, Arg and Lys have 
charge of g = -1-1 e. His has q = +0.98 e, and Glu and Asp 
oi q = —0.76 e. In view of the coarse grained level of the 
model, and to simplify the calculations, we consider that 
each basic residue carries the charge of g = +1 e and each 
acidic residue of q = — 1 e. With such simplifications, the 
total charge of Nl domain is ^ni = +10 e (+15 e and 
-be), for CI it is qci = +9e (+10 e and -le), for NIO 
it is (7N10 = +6e (+8e and — 2e), and the charge of CIO 
is <Zcio = +7e (+8e and — le). The total charge of all 
four terminal domains will be gNC = +32 e. The AA 
model for all four domains is illustrated in Fig. [2] 

C. Polyglycine (PG) model for terminal domains 

When examining the central parts of the residue se- 
quence in the AA model in more detail, one can observe 
a repeating pattern of polyglycine blocks separated by 
one or several H, P or, rarely, basic monomers. One 
can also notice that the acidic residues are mostly lo- 
cated at the beginning of the tails (near the helical IF 
part, represented by planar surface in our SCF model) 
while the basic ones are mostly situated at the end of 
the tails (far from the IF surface). In order to capture 
and emphasize the major specific properties of the ter- 
minal domain structure we have designed a simplified 
"polyglycine" (PG) model for the N and C tails. The 
coarse PG model consists of repeating blocks of four G 
monomers and one H monomer (N tail) or four G and two 
P monomers (C tail) with additional five basic residues 
at the end of each tail. Thus, the structure of the N tail, 
NpG, is Hi[G4Hi]3i(+)5 and the structure of the C tail, 
CpG, is Pi[G4P2]22(+)5- The lengths of the Npg and 
CpG fragments are chosen to be 161 and 138 residues, 
respectively: this is because these numbers are near the 
average of Nl and NIO tail lengths (180 and 146 residues) 
for NpG and, consequently, the average of CI and CIO 
tail lengths (151 and 124 residues) for Cpg. The Npg 
and CpG tail models are also illustrated in Fig. [2j 

D. Modelling parameters 

The short-ranged Flory-Huggins interaction parame- 
ters X between the different types of monomers applied 
for both tail models are the follows. The hydrophobic 
residues, H, strongly repel all the polar ones, so we set 
X = '^k^T for interactions of H with water and ions (Na, 
CI), and X — ^^bT for interactions of H with all the 
polar protein residues (P, +, — ). The interactions of Na 
and CI with water are attractive, x — ~1 ^b?^ to mimic 
the tendency of hydration for the ions. Concerning the 
last residue group, G, we set x = 0-4 for interactions 
of G and H group and x — 0.6 /cbT for those between G 
and all the others residues (water, P, +, — , ions). The 
remaining interactions are set to be athermal (x = 0)- 
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TABLE I: Set of Flory-Huggins interaction parameters x (in 
units of ksT), charges q (in units of e), and relative dielectric 
permittivities er- A dash '-' indicates a zero entry. At (*) the 
surface (s) charge density (in units of a^^) is varied between 
(Ta = — 0.0655 e and as = — 0.071 e for the AA model (at 
(Ts = — 0.0664 e the surface charge is fully balanced by the 
charge of the grafted chains), and between Us = — 0.0405 e 
and o-s = -0.0425 e for the PG model (at = -0.0415 e the 
surface charge is again fully balanced by the charge of the 
grafted chains). 



All the monomer types considered have no affinity to the 
surface, Xs = 0. 

The choice of the interaction parameters for glycine 
is based on the experimental data for solubility of free 
glycine and glycine oligopeptides in water. Experimental 
evidence show that free glycine has rather good solubil- 
ity in water [55] . However, the solubility of oli- 
goglycines is much lower and it reduces with increase of 
the oligopeptide length [13113111]. Lu ei oL [H] mea- 
sured solubilities of glycine and its oligopeptides up to 
hexaglycine at different pH values and the results show 
that the solubility of oligoglycines longer than 3 residues 
strongly decreases with length. Bykov and Asher [17] re- 
ported that oligoglycines longer than 5 residues are nor- 
mally insoluble in water; and Ohnishi et al. (46j stated 
that solubility of polypeptide with glycine linker beyond 
6 is reduced and polyglycine segments longer than 9 
residues form insoluble aggregates. In the current model, 
glycine is present in both forms: as oligomers in the se- 
quence of keratin terminal domains and as free amino 
acid in the NMF composition. Taking into account dual 
hydrophobic-hydrophilic properties of glycine, we antici- 
pate that the interactions of glycine with both hydropho- 
bic and polar residues should neither be strongly repul- 
sive nor attractive. Thus, we set the interactions with 
non-polar residues slightly attractive {x = 0.4 /cb^) and 
with all polar slightly repulsive (x = 0.6 /cbT). Alterna- 
tively, it would be possible to separate glycines into two 
groups: one for free glycine in NMF, and another one for 
glycine blocks in terminal domains, but that is beyond 
the scope of the current simplified model. 

We have considered different values of dielectric per- 
mittivities, Eq, for different species components in our 
calculations. A similar approach has been used by Leer- 



makers et al. for modelling projection domains of neuro- 
filaments [39U43j . The permittivity for water was set to 
Ca — 80, for hydrophobic group H and IF surface — 2, 
for all the polar and charged components (P, +, — , Na, 
CI) £a = 5, and for glycine (G) we set — 4. The 
local dielectric permittivity was calculated according to 
e{z) = Co (-z) J where eo is the permittivity of vac- 

uum and 4>°'{z) is the volume fraction of species type a 
at distance z. The set of all parameters for both A A and 
PG models is given in Table |l) 

The calculations were carried out using the lattice 
spacing of ao = 0.4 nm. There are literature reports of 
lattice spacing ranging between values of 0.3 nm, used for 
modelling of caseins [341 - 138] and 0.6 nm applied for cal- 
culations of terminal domains of NF [55H15] and 3RS tau 
protein [44] . In our model system the main components 
are amino acids and water molecules. The water molecule 
size is about 3.1 A and that of amino acids ranges from 
3.9 A for Gly to 6.1 A for Trp 54J. As glycine, the small- 
est amino acid, is the main component of the sequence of 
the Kl/KlO terminal domains, we used the intermediate 
value of ao = 0.4 nm as lattice size in our calculations. 

In this work we consider terminal domains uniformly 
grafted into two IF cores, which are represented by planar 
surfaces. The grafting density of the domains is calcu- 
lated according to the fact that there are four terminal 
domains (two N and two C) per one dimer length of i = 
50 nm, and one IF core consists of 16 dimers (8 protofila- 
ments) , which gives in total 64 domains per dimer length. 
Taking the IF core diameter of 2R = 7.8 nm [H] and 
the lattice size oq = 0.4 nm, we obtain the grafting den- 
sity a = Mal/{2iTLR) = 0.0083 (in units of ag). For 
the AA model, with the average charge per each tail 
(9n) = (9c) = +8e (as the total charge of the four do- 
mains is (/NC = +32 e), the charge density on the surface 
due to grafted chains would be ctnc — +0.0664 e. As for 
the PG model the charge density due to grafted chains 
would be, correspondingly, ctnc = +0.0415 e. 

The calculation of the IF coiled-coil backbone charge 
is not so obvious due to lack of information about IF 
core organisation. Considering all the charged amino 
acids on the Kl/KlO a-helical parts, we have obtained 
N+ = 45, = 58 for Kl and N+ = 39, iV_ = 60 for 
KIO, which gives the net charge of the Kl/KlO dimer 



— — 34e. We also take into account that 14 salt 



bridges do not change the total charge of the dimer. Car- 
rying out the charge calculations for the IF core with 
dimer length of 50 nm, we consider, as previously, that 
IF backbone comprises 16 dimers in its cross section. 
Thus, we obtained the surface charge density (i. e. charge 
per Cq) of (Js — —0.071 e. This result for the surface 
charge density is appeared to be quite close to the value 
of the surface charge density due to the grafted chains, 
cnc = +0.0664 e (with opposite sign). As explained by 
literature reports [121 [13 113 [56] the IFs are apolar; we 
expect that the overall charge of IF core, dangling termi- 
nal domains, and appropriate counterions should be bal- 
anced. Because the exact IFs organisation is unknown. 
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we can not estimate how many of the accounted amino 
acids on the IF core are in their dissociated form. Keep- 
ing in mind that some parts of the protofilaments and, 
therefore, some of the charged amino acids could be hid- 
den inside the IF core where is no water, we presume that 
the charge density of IF core could be lower then the cal- 
culated value of (7s = —0.071 e. Therefore, as a reference 
(starting) point for our calculations for the AA model we 
consider the surface charge density fully balanced with 
the charge density of the grafted chains, CTs = — 0.0664 e. 
As for the PG model, the balanced value of the charge 
density would be CTs = — 0.0415 e. We also explore a 
range of the surface charge densities around these values 
for both AA and PG terminal models. Salt is represented 
by added Na and Gl ions and the concentration was var- 
ied from as low as Cg = 10~^ M (0s = 3.8 x 10""^) to 
Cs = 0.1 M (0s = 3.8 X lO^''), depending on the system. 

All the calculations were performed using the SGF code 
SFBOX kindly provided by Frans Leermakers. 

III. RESULTS AND DISCUSSION 

A. Overview 

Structural organization of the N and C terminal do- 
mains in one surface and interactions between the two 
IF surfaces will be presented and discussed below. The 
results considered include volume fraction profiles, 0(z), 
of N and G terminals for both models, corresponding 
profiles of the basic residues, 4>+{z), and free energy of 
interactions between the two IF surfaces with the at- 
tached terminal domains, V{D). We start with the case 
of balanced charge densities for the surface and chains, 
Icsl = ctnc, at low ionic strength, Cs = 10~^ M. Then we 
consider the effect of added salt and discuss the options 
when the surface charge in absolute value is higher or 
lower than the charge on the grafted chains. At the end, 
in order to obtain better insights into the properties of 
each type of terminal domains, we consider the interac- 
tions of the surfaces with only one type of the chains (N 
or G) grafted. 

B. IF surfaces and tails at equal absolute charge: 
volume fraction profiles 

The volume fraction profiles, 0(z), show the monomer 
density of the grafted N and G domains at distance z 
from the surface. These distributions provide an estimate 
of how far the grafted chains extend from the surface 
and what is the most probable location of any specified 
monomers. The volume fraction profiles for the whole N 
and G chains are given in Fig. [3] while Fig. |4] shows the 
distributions for only positively charged monomers of N 
and G domains, (ji+i^z). In order to obtain the spatial 
distribution for an unaffected N and G chains, the profiles 
were obtained at large surfaces separation, so that the 



grafted chains do not interact (this corresponds to the 
limit of an isolated IF in solution). 



1. Distribution of N and C tails 

In Fig. [3f^ A) we present the volume fraction profiles for 
the PG model of N and G tails. The profiles for the two 
tails are quite different: the monomer distribution of the 
more hydrophobic N tails is more narrow compared with 
the profile for G tails, with most of the monomers located 
in the first 10 layers from the surface and the maximum 
density at z — Aqq. The extension of the N tails does 
not exceed z — ISag. More polar G tails have lower den- 
sity near the surface and more extended profiles. The 
maximum density is slightly shifted away from the sur- 
face, z = 6ao, and the profiles extend up to z = 20 oq. 
With the simple block-copolymer model for terminal do- 
mains we obtained the two distinct populations of the 
chains: (i) more hydrophobic N tails are collapsed near 
the surface and (ii) more polar G tails are projected far- 
ther into the solution. However, we should notice that 
both types of chains are actually quite compact near the 
surface. With contour lengths of 161 qq and 138 ao the 
chains do not spread out more than 13 oq and 20 oq re- 
spectively. We attribute this behaviour not only to the 
hydrophobic nature of both tails, major component of 
which is glycine, but also to the attraction of the posi- 
tively charged end-monomers to the negatively charged 
surface, causing formation of loops. 

The profiles for more realistic (AA) model for termi- 
nal domains, presented in Fig. |3]^B), show much smaller 
difference between the distributions for N and G tails. In 
general, the behaviour of the all four chains is similar: the 
distributions are quite narrow; most of the monomers are 
located within the first 10 layers from the surface, with 
the maximum density at z = 2-3 qq. The heights of the 
density maxima refiect chain lengths, with the highest 
maximum for the longest Nl tail and lowest one for the 
shortest GIO. Having the contour length of 124-180 oq 
all the chains are in collapsed state and do not protrude 
far into the solution due to their hydrophobicity (^ 50% 
of glycine) and the electrostatic attraction to the sur- 
face. The highest value of volume fraction is obtained 
for Nl, the longest domain (180 aq). With the maximum 
in layer 2, the chains do not extend more than z = 14ao. 
The high monomer density near the surface reflects the 
strong hydrophobic properties of Nl tail — with 24% of 
non-polar and 40% of glycine residues the chains prefer to 
be in compact conformation, reducing contacts with the 
polar solvent. Similar tail extension is observed also for 
NIO except that the maximum density value is lower than 
that for Nl, because NIO chains are shorter and slightly 
less hydrophobic (17% of H monomers and 47% of G). 
The volume fraction profiles for G tails are slightly more 
extended than those for N tails. In particular, the distri- 
bution for GIO tail extends farthest, up to z = 20 oq, and 
the monomer density near the surface is reduced. Even 



6 





5 10 

0-2 -1 ^ AA(N+C) 
/ \ 

// \\PG (N+C) 
;/ \\ 



15 



20 



25 




FIG. 3: Volume fraction profiles for (A) AA model, (B) PG model, (C) comparison of volume fraction for both N tails (Nl+NlO) 
and both C tail (CI and CIO) for the AA model with N and C tails of the PG model, and (D) volume fraction of all the tails 
for the AA and PG models. 



though C tails also consist of about 50% G residues, the 
fraction of non-polar H monomers is much smaller, 9% 
and 2% for CI and CIO, respectively. Being more po- 
lar than N tails, C tails extend a little farther into the 
solution. 

In Fig. and (D) we compare the monomer distri- 
butions for the two models. Fig. [5]JC) shows the profiles 
separately for N and C tails and Fig. [sf^D) compares the 
total profiles for N-l-C tails together. The simplified PG 
model of the N tails gives the density distribution quite 
similar to the combined profile for Nl-l-NlO tails, see 
Fig. |3jC). Even though in the more detailed AA model 
the maximum is slightly closer to the surface and the ex- 
tension of the profile is slightly larger (dashed line), these 
differences are comparatively small. As for the C tails, 
the difference between the profile for the PG model and 
the combined Cl-l-ClO profile for the A A model is more 
pronounced. The general shape of the profiles is similar, 
so is their extension (to z = 20 aq), but the density maxi- 
mum for the PG model is lower and shifted away from the 
surface. That gives the impression that the hydrophilic- 
ity of the C tails in the PG model is somewhat overesti- 
mated; the more accurate AA model predicts that the C 
tails are more hydrophobic. Nevertheless, the relatively 
narrow profiles for the terminal domains coincide with 
the prediction of the compact structure of the tails due 
to formation of the glycine loops [5^. The glycine loops 



hypothesis predict that quasi-repetitive, glycine-rich ter- 
minal domains of epithelial keratins comprise flexible 
and compact glycine loops, where sequences of glycine 
make loops between the stacked non-polar residues. Even 
though SCF method does not allow obtaining such struc- 
tural loops, it predicts compact conformation of the ter- 
minal domains near the surface. Therefore, despite some 
discrepancies in individual profiles for N and C tails, in 
the two models, the combined profiles for all (N-l-C) tails 
are fully consistent with each other, see Fig. [3](D). The 
simple glycine multi-block model for N and C terminal 
domains reasonably well reflects the density distributions 
of terminal domains for Kl/KlO IF. 

2. Distribution of the basic residues of the tails 

The volume fraction profiles for N and C domains pro- 
vide the information about spatial distribution of the 
chains as a whole, while the location of the ends of the 
chains can be obtained from the distribution of the posi- 
tive residues, <j)+{z). The basic residues in the PG model 
are located only at the end of the tails and in the AA 
models they are scattered along the chains with higher 
concentration near the ends. The volume fraction profiles 
for basic residues in N and C tails, (f>^{z), are presented 
in Fig. |4| 
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FIG. 4: Volume fraction of the basic residues on the N and C 
tails for the AA and PG models. 




-0.004 -I 



FIG. 5: Interaction potential V{D) for the AA and PG models 
at the equal absolute charges of surfaces and grafted chains 
and low salt concentration, Cs = 10~^ M. 



For the PG model the distributions of positively 
charged monomers for N and C tails are practically the 
same. The positive monomers for both tails are located 
near the surface, with maximum at the first layer fol- 
lowed by abrupt decrease in the monomer density with 
the distance from the surface. For the distances z > 5 oq 
the fraction of basic monomers becomes very small. That 
result allows us to conclude that the basic residues, and, 
therefore, the end of the chains are located at the sur- 
faces, so the tails form either loops back to the grafting 
surface or bridges with the opposite one. 

As for the AA model, the distributions of positive 
residues for N and C tails differ both from those in the PG 
model and between each other. First, both distributions 
for the AA model are wider, especially for the C tails, 
and second, the difference between the 4>+{z) profiles for 
N and C tails is more noticeable. For the N tails, (t>+{z) 
is similar to that for the PG model, with the maximum 
at the first layer and subsequent decrease of the density 
with distance. At distances z > 10 ap very small fraction 
of basic monomers can be found. The total volume frac- 
tion is higher than that for the PG model because the 



amount of the positively charged monomers is higher. In 
the PG model there are only 5 basic monomers in each 
tail, while for Nl and NIO the numbers of basic monomers 
are 15 and 8, respectively. Taking into account that the 
grafting density of N tails for the PG model is the same 
as the sum of the grafting densities for Nl and NIO, the 
calculated total amount of the positive charges for both 
N tails in the AA model is more than twice higher than 
that for the PG model. That results in about double 
the value of volume fraction of basic monomers for the 
AA model. Positively charged monomers for C tails dis- 
tribute much wider and spreading gradually over ^17 
layers from the surface. The maximal density is again in 
the first layer but its value is more than half than that for 
the N tails, even though the number of positive charges 
for the C tails is not much smaller, 9 and 8 for CI and 
CIO, respectively. 

The density profiles for all the tails show that the maxi- 
mum density for the basic monomers is always at the first 
layer. The fact that the highest concentration of those 
residues is at the surface confirm our hypothesis that 
the charged monomers adsorb onto the surface, so the 
chains form loops and/or bridges between the surfaces. 
Broader volume fraction profiles of basic monomers for 
the AA model possibly results from the different distri- 
bution of the charged monomers along the chains. In the 
more detailed AA model, the basic monomers are not lo- 
cated exactly at the end of the chains, but somehow dis- 
tributed along the whole length of the chains, with higher 
concentration at the ends. Thereby, the more uniformly 
distributed charges in the AA model give a thicker ad- 
sorbed layer while the clustered charges in the PG model 
adsorb fiat on the surface, producing a very thin layer, 
similar to that of highly charged polyelectrolytes. 



C. IF surfaces and grafted tails at equal absolute 
charge: interaction potential profiles 

The interactions between the two surfaces (IF cores) 
covered by grafted N and C terminal domains can be eval- 
uated by calculating free energy of interactions between 
the surfaces at each separation D. The free energy of in- 
teractions A{D) is calculated from the partition function 
under conditions of restricted equilibrium, described by 
Evers et al. |33| . Under such conditions some components 
of the system are free to diffuse from the gap between the 
two surfaces to the bulk solution (e.g. water molecules, 
ions, free amino acids) and the others are restricted to 
stay within the gap (e.g. grafted N and C domains). The 
net interaction potential, V{D), is the difference between 
the free energy value at separation D and its value when 
the surfaces are far apart, V{D) — A{D) ~ A{D^) and it 
is measured in units of kBT/a^. The "far apart" separa- 
tion, Doo, is such that the two surfaces do not interact; in 
our calculations Doo ranges between 150 oq and 1000 ao, 
depending on the model and salt concentration. When 
the interaction potential is negative, V{D) < 0, the two 
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FIG. 6: Interaction potential V{D) for (A) PG model, (B) AA model for various ionic strength indicated on the graphs at 
surface-chains charge balance conditions |o"a| = ctnc- Dashed lines represent data for NMF instead of salt. 



surfaces attract each other, while the positive potential, 
V{D) > 0, implies the repulsive interactions between the 
surfaces. It can be shown that the interaction force be- 
tween the two polymer-covered surfaces can be evaluated 
from the obtained interaction potential [33l [58] . 



1. Low ionic strength 

The interaction potential for surfaces with attached 
N and C terminals (for both models) in conditions of 
charge balance between surface and chains, \as\ = ctnc, 
is presented in Figs. [5] and [6] There is no need for ad- 
ditional counterions to satisfy the charge neutrality con- 
dition. Ideally we would run the SCF calculation in the 
absence of added ions, but we are forced to introduce an 
extremely small concentration of added salt to maintain 
convergence. Still, the extremely low salt concentration 
case captures the experimental set up of two IF surfaces 
immersed in deionised water, where a small salt concen- 
tration cannot be avoided. Such is the case of the Jukura 
experiment [19] where the water extractable materials 
(NMF) from the SC sample were first released and then 
deionised water was added. Other extreme cases, where 
the surface charge is higher or lower than the charge on 
the terminal chains and addition of certain amount of salt 
(counterions) is required to obtain the charge neutrality, 
will be presented and discussed further below. 

Fig. [5] compares the interaction potential for both AA 
and PG models at low salt content, Cg = 10^^ M. For 
both models the interaction potential has a well pro- 
nounced minimum at _D = 17 oq (6.8 nm) for the PG 
model and D = 18ao (7.2 nm) for the AA model, cor- 
responding to net attractive interactions in the system. 
It is interesting that the separations at which the attrac- 
tive minimum occurs {D « 7nm) are in agreement with 
the experimental values for the distance between the two 
IFs, D « 8.2 nm [14]. The attraction between the sur- 
faces at separations D w 15-35 oq (where V{D) < 0) 
occurs due to the well known polyelectrolyte bridging ef- 



fect [30l [33l [Sni [60] and favorable electrostatic conditions 
(ionic strength). In the limit of low surface coverage, 
we believe that positively charged end-monomers are at- 
tracted to the opposite surface forming bridges across 
them. The possibility to be simultaneously attracted 
to more than one surface is more entropically favorable. 
The volume fraction profiles of the charged monomers 
discussed above support this picture as the positively 
charged residues are mostly located near the surface, 
which indicates the possibility of formation either loops 
or bridges (if the surfaces are close enough). At larger 
surface separations, Z? > 35 oq the interaction potential 
approaches zero, indicating that the grafted chains do 
not interact. However, at short separations, D < 15 oq, 
the potential is positive due to strong steric repulsion 
between the chains. 

We should draw attention to the fact that our simple 
PG model for N and C tails, based on the repetitive mo- 
tif of glycine blocks, very well reproduces the result of 
the more complex AA model based on the amino acid se- 
quence. Both characteristics of the system — the volume 
fraction profiles and the interaction potential between the 
IF surfaces — are in a good agreement between the two 
models. We believe that the PG model can be slightly 
improved, for example, by introducing some H residues 
into the C tail model and/or by distributing the charge 
less blockwise along the chain. Despite its simplicity, the 
PG model reflects well the properties of the N and C do- 
mains and, therefore, it probably can be used as a start- 
ing point for more refined (and computational intensive) 
modeling techniques, such as MC, MD or DPD. 



2. High lomc strength and NMF 

The interaction potential for the two models at dif- 
ferent ionic strength is given in Fig. [6] As we already 
discussed, at low salt content, Cg < 0.1 mM, the interac- 
tion potential develops an attractive minimum at short 
separations between the surfaces and levels to zero at 
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H 


13.7% 


P 


25.5% 


G 


10.5% 


+ 


8.4% 




11.2% 


(+)* 


2.8% 


water 


27.9 % 


total 


100.0% 



*Neutraliser. 



TABLE II: Composition of NMF solution. 



longer separations. At higher ionic strength, Cg ~ 1~ 
10 mM, the minimum becomes shallower and a repulsion 
apprears at larger separations. Our two models for termi- 
nal domains give qualitatively similar results but in the 
PG model more added salt is required to destroy the at- 
traction, i.e. the PG model still shows a small attraction 
at Cs = 10 mM, while for the AA model the interactions 
are already repulsive at all separations at = 5 mM. 
That occurs because the attraction to the surface of the 
charged block at the end of the PG chains is stronger 
than that of the AA chains (more uniform charge dis- 
tribution along the chains), so more salt is needed to 
affect the attraction. In every case, at salt concentra- 
tion near physiological conditions, Cg = 0.1 M, the strong 
repulsion between the surfaces with grafted chains is ob- 
tained for both models. High salt content leads to elec- 
trostatic screening, so the grafted chain mediate essen- 
tially a steric repulsion between the two opposing sur- 
faces, much like polymer brushes. Obviously, at even 
higher ionic strength {cg > 0.1 M) the repulsion between 
the surfaces becomes stronger and more short-ranged. 

Calculations at varying levels of salt concentration 
were aimed at mimicking the Jokura experiments with 
normal (healthy) and reduced amounts of natural mois- 
turizing factors (NMF) in SC. The amount of NMF is 
directly responsible for hydration level and elasticity of 
the skin . NMF is made mostly of free amino acids 

derived from the enzymatic degradation filaggrin, as well 
as organic and inorganic salts (19ri26) . As charged amino 
acids and ions are important components of NMF, we 
have used high ionic strength as the first approximation. 
Increasing amounts of added salt tips the IFs interac- 
tion from attractive to repulsive. Experiments show that 
treatment with potassium lactate, could restore the SC 
hydration [20J. 

The next level of increasing complexity in our model 
was to account for the complex mixture of amino acids in 
the suspending matrix between the IFs. In order to cre- 
ate our coarse model of NMF we adopted the amino acid 
composition form Jacobson et al. [25] and then divided 
all amino acids into the same groups (H, P, G, +, — ) used 
in the model for N and C tails. The water content was set 
to 30% as reported in the literature [ini HSJ [51] and 



addition of neutraliser was necessary to ensure charge 
neutrality in the bulk. The detailed composition of this 
NMF -I- water model is given in Table |IT] 

The interaction potentials for two IF surfaces with 
grafted terminals immersed into NMF solution are pre- 
sented in Fig. [6] by dashed lines. The graphs clearly illus- 
trate that the more complex NMF-water mixture leads 
to an even stronger repulsion. We have observed that 
the mixture of free amino acids has stronger effect on the 
interactions between IF surfaces than just adding salt to 
the solvent; the NMF not only provides a strong repul- 
sion between the approaching surfaces but also "pushes" 
the surfaces further away from each other. We believe the 
reason for such a strong repulsion between the surfaces 
rests in the high amount of free charged species (ions 
and amino acids), but what makes this forces more long 
ranged is the presence of free neutral amino acids. Solu- 
tion of only neutral amino acids only slightly decreases 
the attraction between the surfaces and shifts of the at- 
traction minimum to larger separations. 

This result does not support the Jokura et al. find- 
ing that neutral amino acids improve mobility of keratin 
fibers, as well as basic amino acids, but not acidic ones 
[TO] . Our results showed that only charged species in 
solution can affect the attractive intermolecular forces 
between negatively charged IF cores with grafted posi- 
tively charged terminal chains. We should also mention, 
that in our coarse-grained model we could not reveal the 
specific effect of basic amino acids, as the properties of 
positive and negative free amino acids are the same ex- 
cept of the charge and the charge neutrality is required 
in the bulk. In order to examine the effect of specific ions 
a more sophisticated model and/or method is required. 



D. Surface charge higher or lower than the charge 
of the terminals 



Previously we have described the case when charge on 
the surface is fully balanced by charge on the grafted 
chains only. In this section we consider cases when the 
charge of the surface is higher or lower than the charge 
on the grafted tails, so to obtain charge neutrality certain 
amount of counterions is needed. We do so by setting the 
concentration of salt Cg in the bulk, which is in equilib- 
rium with the gap between the decorated IF surfaces. 

In Fig. [7] we present the interaction potential, 
V{D), for the surface charge densities of CTs — 
-0.0655 e, -0.0660 e, -0.0664 e, -0.0670 e, -0.0675 e, 
and -0.0710 e for AA model and cr^ = -0.0405 e, 
-0.0410 e, -0.0415 e, -0.0420 e, and -0.0425 e for the 
PG model. The charge on the grafted chains is kept 
constant at the values of ctnc — 0.0664 e for the AA 
model and (Tnc ~ 0.0415 e for the PG model. Thus, 
for the AA model, surfaces with charge density |(Ts| < 
0.0664 e are "undercharged" (in the specific sense that 
the surface charge density is smaller in absolute value 
than that needed to balance the charge on the grafted 
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FIG. 7: Interaction potential for (A) PG model, (B) AA model for various surface charge indicated on the graphs neutralized 
by grafted tails and ions. The charge of the grafted chains is fixed via grafting density of cr = 0.0083. The amount of salt 
needed for charge neutrality depends on the surface charge and is given in Table [Till 
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TABLE III: Values for surface charge density (cts) and the 
corresponding salt concentration (c^) required for charge neu- 
tralization. The difference between the charge densities of the 
surface and grafted tails is denoted by Aa. 



chains) and, correspondingly, with |(Ts| > 0.0664 e they 
are "overcharged" . For the PG model the threshold val- 
ues of surface charge density for undercharged and over- 
charged surfaces would be, respectively, \as\ < 0.0415 e 
and IcTsI > 0.0415 e. 

For each surface charge we found the values of the salt 
concentration which provide charge balance. The surface 
charge densities with the balancing salt concentrations 



are given in Table III 



In the cases of under- or overcharged surface at low 
salt concentrations, repulsive electrostatic forces domi- 
nate, so the bridging attraction between the covered sur- 
faces could not be seen. At high ionic strength, the re- 
pulsion decreases due to screening. When the surface 
is overcharged, i.e. when the charge on the surface is 
higher in absolute value than the charge of the chains, 
it is possible to find balancing salt concentration, un- 
der which the interaction potential between the surfaces 
would be the same as for the case when the surface charge 
is fully balanced by the charge of the chains only. The 
graphs in Fig. [7] show that at certain amount of added 
salt the interaction potential profiles for |crs| > 0.0664 e 
for the AA model and for \as\ > 0.0415 e for the PG 



model completely overlap. Table [nl] also shows that the 
stronger the charge imbalance (difference between sur- 
face and chains charge), the higher the amount of salt is 
required to neutralize the charge in the system. However, 
when the surface is undercharged, the attractive part is 
reduced and the potential always displays long-ranged 
repulsion, which increases with increasing charge imbal- 
ance. This phenomenology is the consequence of charge 
screening, as the following simplified model calculation 
shows. Consider a plane surface with a negative sur- 
face charge density as, surmounted by a charge cloud at 
charge density ctnc uniformly distributed over a thickness 
H. We solve the linearised Poisson-Boltzmann equation 
for this problem, 



-AttIb<j^c/H {0<z<H) 
{z> H) 



(1) 



where ip is the electrostatic potential (in units of k^T /e), 
Ks is the inverse Debye screening length (k^ — SttIbCs), 
and is the Bjerrum length {Ib ~ 0.72 nm). The bound- 
ary conditions are dip/dz — — 47r/B|crs| at the wall and 
93 — >■ as z — oo, and if should be continuous at z = H 
with a continuous first derivative. This problem can be 
solved analytically. The behaviour of the potential at 
distances z > H from the surface is the relevant piece of 
information, 



CTNC smh{KsH) - \as\KsH 
2Hc, 



X e 



{z>H). (2) 



The prefactor indicates there is a special balance point 
where the potential vanishes completely for z > H. This 
point occurs when IctsI/ctnc — si'nh(KsH)/{KgH). The 
right hand side is an increasing function of UsH, and 
only approaches unity for KsH 0. Thus we see the 
surface has to be overcharged in order to reach the bal- 
ance point and a higher degree of overcharging requires 
a larger value of KgH to compensate, corresponding to 
higher salt, exactly as found above. The reason for this 
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FIG. 8: Interaction potential for N and C terminal domains separately, (A) PG model; (B) AA model. The curves for the cases 
when all the tails are present are also given for comparison (bold lines). 



is that for z > H the surface charge density is screened 
by an factor ^ e^'^'^ relative to the diffuse oppositely- 
charged cloud. 

The interactions between negatively charged surfaces 
covered by positively charged polyelectrolytes were in- 
vestigated experimentally [63, 64] by Monte-Carlo simu- 
lations [64], and theoretically ^65]- The results of these 
studies have shown that the attractive bridging can dom- 
inate only when the charges of the polymers and ions bal- 
ance the charge of the surface. Claesson and Ninham [S3] 
demonstrated that attractive forces between mica sur- 
faces covered by adsorbed chitosan were observed only 
when electrostatic double layer disappeared, i.e. when 
surface charges are exactly balanced by the charged of 
adsorbed polysaccharide. When charge of chitosan, con- 
trolled via variation of solution pH, was higher or lower 
than the charge of the mica surfaces, the electrostatic 
double layer repulsion forces dominate. Dahlgren et al. 
[64] measured the force acting between two mica sur- 
faces covered by MAPTAC polyelectrolyte and also car- 
ried out MC simulations for two surfaces with covered 
by oppositely charged polyelectrolytes. When PE ad- 
sorption was such that the surface charge was balanced 
by the polyelectrolyte, a strong attractive force was ob- 
served at short surface separations. Addition of salt 
to the MAPTAC solution facilitates the increased ad- 
sorption of polyelectrolyte, that leads to a reduced at- 
traction and the appearance of a repulsive double-layer 
force. The authors concluded that the attractive bridging 
mechanism will only dominate when the polyelectrolyte 
adsorption approximately neutralizes the surface charge 
density. Borukhov et al. [65 proposed a theoretical ap- 
proach to explain the behaviour of polyelectrolytes be- 
tween charged surfaces. Their calculations show that at 
low ionic strength the attractive interactions between the 
surfaces take place when polymer adsorption balances 
surface charge. At high ionic strength the surface charge 
is balanced both by polymers and ions and the stronger 
the polymer charge, the more salt is needed to achieve 
the charge neutrality. The authors also considered val- 



ues of adsorbed polymer higher or lower than the equi- 
librium adsorbed amount. When the adsorbed amount 
was lower than the equilibrium one, the attraction was 
weaker. However, when the adsorbed amount was higher 
than the equilibrium one, the results show stronger at- 
traction between the walls and also appearance of strong 
long-ranged repulsion, similar to those shown in Fig. [7] 

In the experiments described by Jokura [IH] loss of elas- 
ticity was observed for SC samples with extracted NMF 
and further hydrated by addition of deionised water. The 
authors suggested that loss of elasticity happens due to 
attractive intermolecular forces between keratin fibers. 
NMF, mainly free amino acids, reduces intermolecular 
forces through nonhelical regions of keratins (N and C 
terminal domains) , so the keratin filaments acquire their 
elasticity. Our modelling results, theoretical consider- 
ation and literature analysis [63^-^65| show that the at- 
tractive interactions between IF at low salt content oc- 
cur only when |(Ts| = ctnc or the IF surfaces are slightly 
overcharged. We conclude that the Jokura experiments 
could take place only at condition that surface charge is 
equal or slightly higher than the charge on the nonhelical 
chains. Thus, the charge of IF cores could not be much 
higher or lower than the charge on the unstructured ter- 
minal domains. 



E. Role of each type of terminal domains 

Why has Nature used two types of unstructured ter- 
minal domains of similar length scale for each keratin 
protein? Does each domain type has a specific function 
and, if so, what is it? Is it necessary to capture the 
specific differences in a model? Different authors have 
taken different approaches. In modelling neurofilament 
projection domains [55H15] . the much shorter globular N 
domains were not included in the study; only the C pro- 
jection domains were considered. On the contrary, in 3RS 
tau protein research ^] the authors focused on the 196 
amino acid long unstructured N domains. Thus, as the 
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role of each domain in keratins is yet unknown, in order 
to generate insights, we decided to take advantage of fast 
computer models and examine the interactions mediated 
by each type of terminal domains separately. 

Fig. |8] shows the interaction potentials for the IF cores 
grafted only with N domains or only C domains. The 
graphs from Fig. [2] summarise the results and compare 
them against the full model calculations. For the calcu- 
lation of only one type of domains, the grafting density 
of the chains was kept the same as before, a = 0.00415; 
this is half the total grafting density for both chain types 
together (a = 0.0083). The charge on the surface was 
then adjusted to neutralize the charge from the chains, 
CTs = -0.0332 e for the AA model and = -0.02075 e 
for the PG model. 

When only N chains are present, the minimum be- 
comes much deeper and is shifted closer to the surface. 
Even with the addition of O.IM of salt this attraction 
minimum is still quite deep (data not shown). Appar- 
ently, the more hydrophobic N tails behave as a "glue", 
holding together the two surfaces. In contrast, the C tails 
behave in the opposite way. The interaction potentials 
for a similar model including C-tails only result in much 
smaller attraction minimum, pushed away from the sur- 
face. The more polar C-tails contribute much less to the 
attraction between the IF surfaces. 

It is tempting to propose that both N and C domains 
play important roles in the structure and interactions of 
skin keratin IFs. The more hydrophobic N chains bring 
about a strong attraction between the IF surfaces while 
the more polar C tails push the surfaces away from each 
other, so that the two types of domains work together to 
keep IFs at the optimal separation. Therefore, we believe 
that it is the combination of both types of the domains 
balances the interactions between the intermediate fila- 
ments. 



IV. CONCLUSIONS 

We have applied the SCF approach to study inter- 
actions of the unstructured N and C terminal domains 
of skin keratin (Kl/KlO) Intermediate Filaments. Posi- 
tively charged N and C domains were grafted onto neg- 
atively charged IF cores, represented by planar surfaces. 
We have considered two models for N and C tails, char- 
acterized by a different level of detail: the coarse block- 
copolymer PG model and the more detailed AA model, 
which is based on the amino acid sequence. In spite of 
the apparent simplicity of the PG model, it qualitatively 
captures most of the effects observed for the more com- 
plex AA model. We have presented monomer density 
profiles for the N and C tails and, separately, profiles for 
their basic residues only. We have compared and dis- 
cussed interaction potential profiles for IF surfaces with 
attached tails at various surface charge densities, ionic 
strengths, and for the solution of free amino acids repre- 
senting NMF. We have also attempted to clarify the role 



of each type of terminal domains considering N and C 
chains separately. Our main findings are summarized as 
follows. 

(A) Volume fraction profiles for N and C domains show 
that the monomers of both types of the chains are mostly 
concentrated near the surface, so the chain extension does 
not exceed r « 20ao = 8nm (Fig. [3]). The basic residues 
of the terminal domains, which are located near the end 
of the chains, have the highest density at the (oppositely 
charged) surface (Fig. |4]). These results indicate that the 
chains form either loops or bridges with another surface. 
Such bridges lead to attractive interactions between the 
two IF surfaces at short separations. N tails are more 
hydrophobic and the profiles for N tails are more narrow 
compared to those for C tails and extend for no more than 
r « 15 oq from the surface. The interaction potential for 
surfaces covered by N domain type only reveals that the 
attractive interactions between the surfaces are stronger 
than those when both types of the domains considered 
together and appear at shorter separations (Fig. [s]). So 
we conclude that N tails work as the "glue" between IF 
surfaces. C tails are slightly more polar than N tails and 
extend shghtly further into the solution (r sa 20ao). The 
interactions between surfaces with only C tails grafted 
show much weaker and more long-ranged attraction. So, 
we propose that C chains are "responsible" for keeping a 
certain distance between IF. Hence, each type of the ter- 
minal domains has its specific role and their combination 
retain IF at certain distance. 

(B) When the charge of the IF surface is neutralized by 
the charge on the grafted chains and the ionic strength is 
low, IFs experience attractive force between each other at 
surface separations D w 15-35 ao (6-14 nm) due to bridg- 
ing effect of grafted terminal domains. This attraction 
becomes weaker and turns into repulsion with increase of 
ionic strength as a result of electrostatic screening. The 
repulsion become stronger and longer ranged when sim- 
ple aqueous electrolyte solution between the IF surfaces is 
replaced by a complex "broth" of amino acids — a coarse 
grained representation of NMF in 30% water. However, 
we can not confirm experimental observations of Jokura 
et al. jl9| that neutral amino acids alone produce a simi- 
lar effect. We have found that charged small species such 
as ions or charged amino acids are necessary components 
of NMF and their role is to decrease electrostatic forces 
between IF. The effect of salt when ions differ not only 
by their charge but also by size should be investigated 
using more complicated model for salt molecules or/and 
by other simulation methods. 

(C) At low ionic strength the attraction between the 
IF surfaces can be obtained only when the charge on the 
surface is fully compensated by the charge on the chains 
and ions. That occurs only when the surface charge is 
equal or slightly higher than the charge of the grafted ter- 
minal domains. Therefore, we propose that: (i) negative 
charge of the IF helical part is equal in absolute value or 
slightly higher than the positive charge of the IF termi- 
nal domains; (ii) the function of NMF is to prevent the 
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attractive forces between protruding terminal domains 
and IF helical cores. When NMF are removed or their 
amount is highly reduced these attractive forces "glue" 
keratin Intermediate Filaments and reduce the elasticity 
of the corneocytes. 
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